Mem. S.A.It. Vol. 0, 1 
©SAIt 2004 



Memorie deiia 




Sh2-188: a model for a speedy PN 

C.J. Warring 1 , TJ. O'Brien 1 , A.A. Zijlstra 2 and J.E. Drew 3 

' Jodrell Bank Observatory, School of Physics and Astronomy, The University 
of Manchester, Macclesfield, Cheshire, SK11 9DL, United Kingdom, e-mail: 
cwareing@jb . man . ac . uk 

2 School of Physics and Astronomy, The University of Manchester, Oxford Road, 
Manchester, Ml 3 9PL, United Kingdom 

3 Imperial College of Science, Technology and Medicine, Blackett Laboratory, Prince 
Consort Road, London, SW7 2BW, United Kingdom 



Abstract. Sh2-188 is thought to be an ancient planetary nebula in the galactic disk. It 
appears to be one-sided with recent observations revealing structure behind the filamentary 
limb. We postulate that Sh2-188 is interacting with the ISM and simulate it in terms of a 
"triple-wind" model comprising of the usual "fast" and "slow" interacting stellar winds plus 
the wind due to motion through the ISM. We have run simulations at various velocities of the 
central star relative to the ISM and find that a high velocity of 125 kms~' best approximates 
the observed structure. We also suggest that Sh2-188 is younger than previously thought 
and that much of the mass lost on the AGB has been swept downstream. 
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1. Introduction 



The planetary nebula Sh2-188 l lSharplessI 
119591) has a one-sided, semi-circular 
and filamentary appearance. It is lo- 
cated in the Galactic plane at I = 128° 



4°. Spectroscopic work (tParke 



lLozinskava and Esipo vl U97lt Qohnso 



1964; 



1975) 



revealed an extremely high NII/Ha line-ratio, 
thought to be due to an overabundance 

of nitrogen. Furth er spectroscopic work 

(Rosa do and Kwitterll 19821) found abundances 
simila r to Peimbert Ty p e I planetary nebulae 
(PN) JPeimbertl 119811) . iKwitter etal.l i 19881) 
presented an identification of a candidate 
central star of the nebula using geometric, 
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colour and apparent magnitude methods. The 
candidate star is not at the geometric centre of 
the optical filaments. 

The distance to the candidate star has been 
estimated at d = 965 + l°n° pc by non-LTE mod- 
elling Napiwotzki ( 1999, 2001) which also es- 
timated its temperature to be 102000 + 32000 
K. Consistently, this nebula is found to be 
an evolved object. iTweedv and Kwitten ( 1996) 
published an atlas of ancient PNe including 
Sh2-188. 

Figure ^ shows an image of Sh2-188 
created by a mosaic of Ha observations 
taken as part of the Isaac Newton Telescope 
Photometric Ha Survey of the N orthern 
Galactic Plane (IPHAS) jDrew et alJ 120051) . 
The bright portion of the nebula in the south- 
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Fig. 1. An image of Sh2-188 created by a mosaic of 
images taken as part of the IPHAS survey in 2003 
(see text). North up and East to the left. The image 
is 9.4 x 8.5 arcminutes. 



east describes a semi-circular arc with a diam- 
eter of 10 arcminutes. Stru cture first noted by 
iTweedv and Kwitteil i 19961) to the northwest of 
the nebula h as now been reveal ed by the ob- 
servations of iDrew et all ( 120051) to be a faint 
circular closure of the bright nebula with con- 
necting tails stretching behind. Figure|2]shows 
this structure. From the front of the filamentary 
structure to the connecting tails the nebula is 15 
arcminutes long. The tail seems to be stretch- 
ing away from the nebula in opposition to an 
inferred direction of motion from the bright arc 
which could explain the one-sided brightening 
as an interaction with the ISM. 

The accepted theory of PN formatio n is the 
intera cting stellar winds model (ISW) JKwokl 
1982) where a fast wind (~ 10 3 kms -1 ) from 
the hot central star of a PN blows into the 
slow wind (~ 10 kms 1 ) produced during the 
preceding AGB phase resulting in the inner 
regions of the slow wind being compressed 
into a dense shell and ionised by the energetic 
UV radiation of the central star. The familiar 
ring-like appearance of PNe is then observed. 
Interaction with the ISM was thought to be 
somet hing that was only obs erved at late times 
( ITweedv and Kwitteil fl9%l) . IVillaveil J200l 
(hereafter referred to as VGM) reopened the 
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Fig. 2. Same as FigureQbut with scaling chosen to 
reveal further structure. North up and East to the left. 
The image shows a faint ring completing the bright 
semi-circular arc and tails stretching away and con- 
necting in the north-east. The image is 18.8 x 17.0 
arcminutes. 



PN-ISM interaction debate by performing sim- 
ulations which follow the full evolution of the 
AGB phase followed by the PN phase with a 
conservative proper motion of the central star. 
Their results showed that the AGB wind struc- 
ture is strongly influenced by the ISM interac- 
tion and thus changes the appearance of the PN 
at early as well as late times. 

Sh2-188 appears to be one of the most ex- 
treme cases of PN-ISM interaction, made par- 
ticularly rare by its filamentary nature. We have 
developed a "triple-wind" model adding a third 
wind due to movement through the ISM into 
the ISW model. We are investigating whether 
this triple-wind model can reproduce the neb- 
ular shape of Sh2-188 using a new code devel - 
oped by Wareing and O'Brien ( Wareina 2005). 

2. The numerical method 

The numerical scheme used to solve the hy- 
drodynamic equations use s the s econd -order 
Godunov scheme due to iFalld dl99ll) . The 
scheme is in 3D cartesian coordinates, fully 
parallel and includes the effect of radiative 
cooling. A small amount of viscosity is re- 
quired in order to control the effects of the 
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Quirk instability llOuirkl fl992). The correction 
takes the form of an addition to the solution 
of the Riemann problem. This additional flux 
is small in regions of smooth flow and affects 
only regions containing discontinuities. 

The numerical domain consists of a grid 
with 200 cells in each direction and uniform 
grid spacing. The central mass-losing star is 
placed at particular coordinates (xo,yo,Zo) an d 
simulations are performed in the frame of ref- 
erence of the star. Mass-loss is effected by 
means of artificially setting the values of den- 
sity, momentum and energy density in the cells 
in a volume-weighted spherical region of ra- 
dius 5 1 cells centred on the star. The condi- 
tions within this region are reset at the begin- 
ning of every computational timestep so as to 
keep driving the wind. Movement through the 
ISM is set to be parallel to the x-axis and ma- 
terial flows in at the (x = 1, y, z) boundary face 
with a positive v x . All other boundaries have 
free-flow conditions allowing material to flow 
out of the domain freely. 

3. The input parameters 

The ISW model of the PN-ISM interaction fol- 
lows the evolution of the AGB and post-AGB 
phases from the beginning of the AGB phase 
through the post-AGB or PN phase. In the 
triple-wind model, this evolution is combined 
with a constant movement through the ISM. 
The evolution is followed in the frame of the 
central star emitting the stellar winds as the 
ISM flows through the computational grid with 
particular velocity Vism- 

Mass-loss via a stellar wind has been 
modelled as a spherically symmetric constant 
mass-loss rate M with constant velocity v. The 
density in the stellar winds is set according to 
p oc r a . In the AGB and post-AGB winds a — 2 
and the constant of proportionality is M/4 n v. 
In the ISM wind a — and the density and 
velocity are constant. 

Typical estimates from the literature are 
used for the wind inputs. The fast stellar wind 
parameters are denoted by a subscript fw: 
M fw = 5 x 10- 8 M yr-\ v fw = 1000 kms^fc 
7f w = 5 x 10 4 K. The slow stellar wind parame- 
ters are denoted by a subscript sw: M sw = 10~ 6 



M B yr-\ v sw = 10 kms^& T sw = 10 4 K. The 
switch between the AGB wind and the post- 
AGB wind is instantaneous and occurs after 
10 years of AGB evolution, which is short 
but allows a stable structure to develop. The 
detailed properties of these winds, whilst be- 
coming more certain, are still too uncertain to 
justify more detailed temporal modelling. We 
are aiming at developing a simple model to in- 
vestigate the ISM interaction. The ISM itself 
is assumed to be homogeneous with character- 
istics of warm neutral medium, the main con- 
stituent of the ob served ISM: 7 \sm = 2500 K 
& Pism — 1 cirT 3 (Burton 1988). The gas pres- 
sure in all three winds is calculated ignoring 
magnetic pressure via P — p k T /m a where k 
is Boltzmann's constant and m. d is the average 
mass of a particle. 

4. Results 

Various values of Vism have been considered 
from kms~' to 175 kms~' in steps of 25 
kms -1 . The evolution o f these models can be 
found in Wareing ( 2005). We present here only 
the case of 125 kms -1 which we find best re- 
produces the general morphology of Sh2-188. 

Figure |3 shows a density slice through the 
computational grid (x,y = 100, z) 2000 years 
after the onset of the fast wind. This snapshot 
has been selected for its resemblance to Sh2- 
188. The gas structure shows a strong bow- 
shock in the direction of motion, a faint shell 
of compressed AGB material completing the 
circle and connecting tails downstream of the 
main nebula made up purely of AGB material. 
Considering the evolution of the gas up to this 
point, the AGB wind has formed a bow-shock 
interaction upstream of the nebula. The com- 
pression and deformation has formed a high 
density bow-shock which reached a stable state 
between the material supplied by the AGB 
wind and ram-pressure stripping of AGB ma- 
terial downstream after 60, 000 of the 100, 000 
years of AGB evolution. The material we see 
in the tail is purely AGB material and thus far 
older than the nebula. During the PN phase the 
wind decreases in density and increases in ve- 
locity resulting in an adiabatic shock that forms 
the bright shell of the nebula. Initially this 
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Fig. 3. A slice through the computational grid at 
(x,y = 100, z) showing the density of the gas 2000 
years after the onset of the fast post-AGB wind. The 
scaling is logarithmic. 



shock is within the region of undisturbed AGB 
wind inside the inner shock and the nebula ap- 
pears circular. After 2000 years, the fast adi- 
abatic shock has expanded enough to interact 
with the bow-shock and strong asymmetries 
appear, similar to the appearance of Sh2-188. 
After this time, the nebula departs from a cir- 
cular appearance as the adiabatic shock moves 
downstream through the undisturbed AGB ma- 
terial. 

5. Discussion 

We have found that the triple-wind model can 
reproduce the appearance of the nebula Sh2- 
188. If we accept the model, the nebula would 
appear to be moving at a high speed of 125 
kms through the ISM. The typical velocity 
of a central star which results in a PN has 
been shown to be in the area of 35 - 50 kms 1 
llBinnev and Merrifieldll99l . Thick disc stars 
resulting in a PN have a highe r average veloc- 
ity of 50 - 75 kms 1 (Binnev and Merrifield 
1998). The central star of Sh2-188 is therefore 
a rare object for it to be moving this quickly. 
Having investigated other proper motions and 
the resulting nebula shapes, a velocity of less 



than 75 kms" does not result in connecting 
tails behind the nebula. At 100 kms 1 fluid 
instabilities seem to destroy the smooth tail 
structure we observe and thus 125 kms -1 is in 
fact the lowest speed at which we see a nebula 
resembling Sh2-188. 

The general idea of a PN moving 
thr ough the ISM was discusse d in depth 
by iD gani and Sokeil ill 994 Il998l) . Rayleigh- 
Taylor instabilities are thought to be able to 
fragment a bow-shock in the direction of mo- 
tion. Magnetic fields are shown to be able to 
alter the fragmentation. In the case of Sh2- 
188, it is thought to be in the galactic plane 
and thus t he observed fragmentation can be in- 
terpreted (Daani and Soket][l998|) as being an 
effect of flu id and magnetic field in stabilities. 
Importantly, Dg ani and Sokerl ill 9981) conclude 
that fragmentation of a substantial part of the 
halo requires a minimum proper motion of 100 
kms -1 . This limiting speed would predict that 
if Sh2-188 has a proper motion of 125 kms -1 
then the bow-shock would fragment. Whilst we 
do not see this in the simulations, we do ob- 
serve this in the images e.g. Figure^ The lack 
of fragmentation in our simulations may be the 
result of a lack of resolution in the computa- 
tional grid. 

VGM pointed out that the problem in avail- 
able models of PN-ISM interaction is that the 
interaction had only been studied by consid- 
ering the relative movement when the nebular 
shell was already formed. The 2D simulations 
they performed for typical PN conditions (ISM 
density: 0.1 cm" 3 ; proper motion: 20 kms" 1 ) 
followed the full AGB and post-AGB evolution 
of the stellar wind. They found that high veloc- 
ities, high ISM densities and magnetic fields, 
some or all of which had previously been re- 
quired to produce observed ISM interactions, 
were no longer necessary to explain the ob- 
served asymmetries in the external shells of 
PN. Further, the stripping of mass downstream 
during the AGB phase provides a possible so- 
lution to the problem of missing mass in PN 
whereby only a small fraction of the mass emit- 
ted during the AGB phase is inferred to be 
present during the post-AGB phase. The shape 
and appearance of our simulated PN is consid- 
erably affected by the evolution of the stellar 
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wind during the AGB phase. We have shown 
that at high proper motions, the ISM interac- 
tion is apparent within a few thousand years of 
the onset of the fast wind and is strongly influ- 
enced by the preceding AGB evolution. 

There is now little doubt that interaction 
with the ISM considerably alters the amount 
of mass within the circumstellar envelope dur- 
ing the AGB and post-AGB phases. Much of 
the mass emitted from the central star is down- 
stream and in fact forms some of the structures 
we observe. In our simulation we find that hav- 
ing injected 0.102 M of material into the neb- 
ula region, only 0.005 M remains and 0.096 
M has been swept downstream. It is impor- 
tant to note that these results are for only one 
set of stellar wind and ISM parameters. Tuning 
the parameters can produce larger nebula radii 
which take longer to reach a point of stability 
and thus more material is in the nebula region. 
This also makes it difficult to place a constraint 
on the distance of Sh2-188. 

The faint emission completing the ob- 
served ring is particularly important. In the 
simulations this structure has a transitory na- 
ture and moves quickly away from the central 
star. This faint ring may be a very good tracer 
for the age of the nebula. The presence of the 
ring close to the star and still appearing to com- 
plete the circle of the nebula leads us to suggest 
that Sh2-188 is of an intermediate age, not an 
ancient object. 

6. Conclusions 

We have simulated the formation of a PN 
where the central mass-losing star moves 
through the ISM and have found that we can 
reproduce the general shape of the nebula Sh2- 
188. Following the AGB evolution has been 
the only way to recreate the whole structure, in 
particular the tail behind the nebula compris- 
ing of AGB material alone. ISM interaction in 
PN is observable from an early to intermediate 
age; Sh2-188 is most likely of an intermediate 
age, indicated by the position of the transitory 
closed ring which rapidly evolves downstream 
at the start of the PN phase. The ISM interac- 
tion can explain the problem of missing mass 
in PN and simulations can quantify this effect. 



The next generation of telescopes, particularly 
ALMA, will reveal cool dust structure in the 
universe and shed light on the AGB material 
swept downstream of moving PN. 
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